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Abstract 


The quasi-elastic scattering of two nuclei is considered in the high- 
energy optical model. Energy loss and momentum transfer spectra 
for projectile ions are evaluated in terms of an inelastic multiple- 
scattering series corresponding to multiple knockout of target nucle- 
ons. The leading-order correction to the coherent projectile approx- 
imation is evaluated. Calculations with uncorrelated wave functions 
are compared with experimental results. 

Introduction 

The assessment of radiation risk to astronauts from cosmic radiation is currently an area 
of active investigation. Predictions of biological damage will ultimately require a knowledge 
of the particle fluence spectra at the endpoint of interest. In turn, these particle fluence 
spectra are determined from charged-particle transport codes that must contain a description 
of all important physical processes that occur as the incident ions and subsequent generation 
fragment nuclei pass through natural and protective radiation shielding. A theoretical model 
for the prediction of fragmentation cross sections is extremely useful, as it cannot be expected 
that enough experiments will be performed for all the collision pairs and energies of interest in 
cosmic ray studies. Experimental data are most often in the form of inclusive measurements 
where a single reaction product is typically detected in a many-particle final state. Although 
several mechanisms may lead to the single product, models must be compared with the inclusive 
measurements for validation. 

The inelastic collision of two nuclei at intermediate or high energies is often described as 
a two-step process. The first step includes multiple scatterings between projectile and target 
nucleons leading to the knockout of nucleons and clusters, the production of particles, and the 
deposition of energy. The second step involves the cascade of initially struck particles within their 
host nuclei and the de-excitation of the nuclear systems, which may proceed through particle 
emission. Recently, we have considered inclusive heavy-ion scattering using the high-energy 
optical model (refs. 1 arid 2). The multiple scatterings between projectile and target nucleons 
can be divided into elastic and inelastic collision terms corresponding to a distortion effect and 
the knockout of nucleons, respectively. Calculations (ref. 2) with the independent particle model 
(IPM) show that even for large collision pairs, the number of inelastic collisions that occur is 
quite small, usually two to four. Although correlation effects may be important corrections 
to the IPM, especially for cluster knockout, the rapid convergence in the number of inelastic 
collisions favors a “doorway” picture of heavy-ion fragmentation. In the doorway picture the 
first step involves only a small number of knocked-out. nucleons, with the subsequent motion 
of these particles and nuclear de-excitation leading to the large number of final fragmentation 
channels observed in experiments. 

In our previous work (refs. 1 and 2) the cross-section distribution in total momentum transfer 
was considered. Herein, we extend this work in order to evaluate the energy loss cross section 
of the projectile in inclusive scattering. The relationship of the projectile energy loss to the 
target response function and excitation energy is considered and the effects of multiple inelastic 
scatterings are treated. The correction to the coherent projectile assumption (ref. 1) is evaluated 
to leading order. Previous calculations in high-energy formulations (refs. 3 to 5) have considered 
only elementary projectiles and usually assume a zero-range two-body interaction (ref. 4) or 
a factorization approximation (ref. 3). In this work only forward-peaked wave functions are 
assumed in the evaluation of higher order terms. The model presented herein is immediately 
applicable to the development of nuclear transport codes for bulk shielding materials, as is 


illustrated with calculations of energy loss distributions for protons, 4 He and I6 0 in common 
shielding materials. The methods developed herein are expected to lead to an improved 
description of the first step of heavy-ion fragmentation. 


Multiple Inelastic Collision Series 

In the eikonal coupled channels (ECC) model (refs. 6 to 8) the matrix of scattering amplitudes 
for all possible projectile-target transitions is given by 

7 (q) = ^J Sb e* q b [e*( b ) - l] (1) 

where barred quantities represent matrices, b is the impact parameter vector, q the momentum 
transfer vector, and k the projectile-target relative wave number. In equation (1), Z is an order- 
ing operator for the z-coordinate, which is necessary only when noncommuting two-body inter- 
actions are considered. The eikonal phase elements are defined by matrix elements of arbitrary 
projectile-target states of the following operator (for commuting two-body interactions): 

*< b ) - s^E/ (o) w 


where a and j label the projectile and target constituents, respectively, s is the projection of 
the internal nuclear coordinate onto the impact parameter plane, is the nucleon-nucleon 
scattering amplitude, and k^ is the nucleon-nucleon relative wave number. 

In treating inelastic scattering we assume the off-diagonal terms in y, denoted by xo> are 
small compared with the diagonal ones, \ D , and then expand / in powers of \Q : 

7 (q) = J dhe^e^ib) £ |M!!j (3) 
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We also will make the assumption that the diagonal terms are all represented by the ground- 
state elastic phase x- Using equation (3) we sum over target final states X (continuum) to find 
the inclusive angular distribution for the projectile when its mass remains unchanged as 
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Equation (4) allows only for a study of the momentum transfer spectra of the projectile. In 
considering the projectile energy loss, energy conservation must be treated. With continuum 
states used for the target final state, energy conservation leads to 
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where Epi is the energy of the projectile in the final state, w is the projectile energy loss, and 
we define 
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where ky is the wave number vector of a knocked-out target nucleon. 

The inelastic collision series of equation (5) is expected to converge fairly rapidly, and in the 
next section we consider evaluating this series for an uncorrelated target wave function and using 
plane-wave states for knocked-out nucleons. First, we briefly consider the lowest order term and 
its relationship to the target response function. The first inelastic term is 

Wi (b, b', w) = J d 2 q d 2 q e iq b e~ tq ' h> F (q) F (-q') / NN (q) /jJ N (q') 

X / (^2 6 ( E f ~ E i) < °t\ e * q Sj k >< k | E, e Q &j \°T > ( 8 ) 


If we neglect target recoil we can write 
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and then 
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where the target response function is 
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Only the f = j terms contribute for nonzero w , and we neglect the j / j f terms since bound 
states dominate for small w . The plane-wave impulse approximation (PWIA) assumes that 
X ~ 0 and that only the m — 1 term contributes, such that 
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We also note that in nuclear matter, translation invariance implies (ref. 9) 

(q, q', w) = i? (q, q', w) <5 (q - q') (13) 

which leads to some simplification in distorted wave calculations. Bertsch and Esbensen (refs. 10 
to 12) have introduced a surface response model in which the PWIA is assumed; however, eikonal 
waves are introduced into the response operator and provide localization to the nuclear surface. 
This approach works quite well when multiple inelastic collisions do not contribute. For large 
momentum transfer and larger projectiles these higher order terms may become important and 
are treated next. 

Collision Terms 

In evaluating the collision terms W m) we assume an uncorrelated wave function for the 
target and plane waves for continuum states. The inclusion of final-state interactions occurs in 
the transition form factors of the target discussed below. The projectile motion is treated in the 
coherent approximation following reference 1 and the leading-order correction is considered. 

The first collision term is written as 

Wl (b,b» = - f^ A \ 2 [ d 2 qd 2 q'e^ b W F (q) F (q') 

X /NN (q) 4n fa') J ^2 6 ( W “ E k) Go T k (q) gI Qt (q') (14) 

It is helpful to change variables as follows: 


a = \ (q + q') 

(15) 

P = q - q 

(16) 

X = s — s' 

(17) 

y = \{ s + s ') 

(18) 

XI 

1 

II 

(19) 

S = i(b + b') 

(20) 


such that 


w > <r ’ s ' -> - / rf2 “ ^ ^ (“ + f ) ^ (° - f ) 


(27T% N ) 
x R\ (a,/3, w) 


where we have defined 
and 


A{q) = F (q) / NN (q) 

Ri (a, (3,w) — J 6 ( w ~ E k) Go T k ^ G^ (a~P 


( 21 ) 

( 22 ) 

(23) 


4 



Following Krimm et al. (ref. 3), we can formally treat the delta function in equation (23) by 
introducing a Fourier transform pair: 
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For a nonrelativistic nucleon we have 
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where eg l is the binding energy. Equation (26) then becomes, with the assumption of plane 
waves for the target final state in Gq t ^, 
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where $ is the single-particle wave function of the target ground state. After two integrations 
we find from equations (28) and (24) that 
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The second collision term is found in a similar fashion to be 
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With similar coordinate changes as those described above, the mth order collision term is 

W ro (R,S,ti;) = 



2 


where 


X R m (qi , C*2, ■ • • , atm, 01,02,---, 0m, w) 


r&/nr=i 

x $ (yj + %) <1>f (y> - if-)] 


(32) 


ftm («i , «2, - - - , a m ,/3i , 02, - ■ • , £ m , w) = S 




2m. v (® - ffl m )X)7=l I J 


I Cm— 1}/2 


x Jm — 1 
0 


(\/ 2m,v (“ “ C «J Ej=i x j) 


(ui>e Bm ) 

(“' < c O 


(33) 


We next look for a simplification of the terms for m > 1. 
The mth order Bessel function is given by (ref. 13) 
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and 
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Equation (35) is expected to be a useful approximation since the wave functions are peaked at 
Xj = 0. Also, since e Bl < e # 2 < ... < e Bm and successive terms in the collision series dominate 
as w increases, should not be too large in the region of interest. We then have for the energy 
loss spectra (eq. (5)) 
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The coherent approximation assumes the projectile remains in the ground state throughout 
the scattering. The leading-order correction to the coherent terms occurs in Vi '2 and corresponds 
to the following replacement (from eq. (55) of ref. 1): 



which physically allows the projectile to dissociate in the intermediate state. Further modifica- 
tions, which are not included herein, are necessary when correlation effects are treated. Next we 
consider model inputs and application of the above formalism. 


Calculations 

We next discuss physical inputs necessary for evaluation of the cross sections of equations (39) 
and (40). We employ a two-body amplitude of the form 

/NN (q) = gig+OM e -B, 2 /2 (42) 

47 r 

where the spin-isospin averaged energy- dependent parameters are the two-body cross section a , 
the ratio of the real to imaginary parts of the forward amplitude p, and the diffractive slope 
parameter B. 
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For nuclei with mass number A < 16 we use harmonic oscillator shell model wave functions. 
For s-shell nucleons (ref. 14), 


i \ 3/ 4 

(r) = ( — 1 P -r 2 /^ 


na 


) 


and for p-shell nucleons 
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where a is related to the nuclear radius R by a ~ R 2 and the internal coordinate is expressed in 
terms of spherical coordinates as 
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The s-shell and p-shell probabilities are given by 
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with C s = 1 and C p = 0 for A < 4. 
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The projectile form factors are evaluated from equations (43) and (44) as 

^(q) = j ^C,, P i^(r) | 2 Sr 
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The functions W\ (eq. (21)) and R\ (eq. (29)) are now found after some effort. 

Results and Discussion 

Experiments with 800-MeV protons (ref. 15) were performed with several targets to study 
the quasi-elastic peak and extend to regions of energy loss corresponding to pion production. 
Herein we consider only the quasi-elastic and dip regions of the data. In figures 1 to 6 we show 
comparisons of calculations with experimental data for ^Li and ^ 2 C at several scattering angles. 
The dotted line is for the single- knockout term, the dashed line is for the single- and double- 
knockout terms, and the solid line includes all contributions up to m — 4. The fourth-order term 
makes only a small contribution for both targets. Because of our neglect of target recoil and any 
target excitation energy, the position of the peak of the calculated values is shifted to slightly 
lower energies than those of the experimental values. Following reference 4 we corrected this 
shortcoming by performing calculations at w -> w - e B{ . The position of the quasi-elastic peak 
is then well reproduced by the calculations. In the dip region, pion production channels prevent 
a direct comparison with the data at this energy. The spin dependence of the nucleon- nucleon 
amplitudes also affects the results in this region. 

In figures 7 to 10 we show comparisons of calculations at various angles for 4 He- 4 He scattering 
at 1.05 GeV/amu with data from reference 16. The solid line shows the m = 1 to 4 terms with the 
m — 2 term correction for incoherent projectile motion of equation (41) included. The dashed 
line neglects this correction, the dash-dot line is for the m = 1 and 2 terms, and the dotted line 
is for just the m = 1 term. The coherent projectile assumption decreases the contribution from 
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the second collision term. The comparisons with experimental values in figures 9 and 10 are at 
scattering angles corresponding to momentum transfers where the Gaussian wave functions used 
are known to have insufficient strength. In figure 11 we show the angle-integrated cross section 
for a-ct scattering. The dashed line is the single knockout and the solid line is the complete 
series. Multiple nucleon knockout represents a large correction for large energy loss. 

In figures 12 and 13 we show inclusive 4 He scattering on ie O at 1 GeV/amu for scattering 
angles of 1° and 4°, respectively. The dashed line is for the first inelastic collision term only, 
the dash-dot line is for the first and second terms, the dotted line is for the first to third terms, 
and the solid line is for the sum of the first to fourth inelastic collisions. The higher order terms 
are more important here than they are for the case of proton projectiles. (See figs. 1 to 6.) In 
comparing figures 12 and 13 we note that the forward peak in the cross section is an indication 
that the projectile is unlikely to receive both a large energy loss and a momentum transfer 
without suffering a change in mass. In figures 14 and 15 we show a similar comparison for ^’O 
scattering on 1 f, ( ) at 1 GeV/amu for scattering angles of 0.5° and 1.0°. 

Conclusions 

The high-energy optical model is used to describe energy loss spectra of projectile nuclei in 
high-energy collisions. An inelastic multiple-scattering series is found for inclusive projectile 
scattering that corresponds to the knockout of target particles. Preliminary calculations are 
presented for proton, 4 He, and 1 11 0 projectiles with an approximation to the higher order 
(>2) inelastic collision terms. Improvements in the model should be the inclusion of final- 
state interactions of knocked-out target nucleons and the addition of pion production into the 
calculations of the inelastic spectra. Calculations will also be improved by considering spin- 
dependent two-body amplitudes and the use of response functions that treat low-lying collective 
behavior of the target. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
April 9, 1992 
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Symbols 


A P 

Ap 

b 

c 

E f 

Ei 

Ev 

Ept 

F 

f NN 
G 

J 

k 

k 

fc N N 
mjv 
|O p > 
I Ot > 

P 

q 

R 

R P 

Rp 

r 

S 

s 

X 

y 

a 

cx 

(3 

6 


mass number of projectile nucleus 
mass number of target nucleus 
impact parameter vector 
speed of light 
total final energy 
total initial energy 
energy of outgoing nucleon 
energy of projectile in final state 
projectile one-particle form factor 
nucleon-nucleon scattering amplitude 
target transition form factor 
target constituent index 
projectile-target relative wave number 
wave number vector 
nucleon-nucleon relative wave number 
nucleon mass 

projectile initial state vector 
target initial state vector 
momentum 

momentum transfer vector 
= b-b' 

projectile matter radius 

target matter radius 

internal nuclear coordinate vector 

= 2 (b + b') 

projection of internal coordinate onto impact parameter plane 
= s — s ; 

= 5 ( s + s') 

projectile constituent index 

= 5 (q + q') 

= q-q' 

Dirac delta 
binding energy 
cross section 
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$ single-particle wave function of target ground state 

X ground-state elastic eikonal phase 

X ground-state eikonal phase operator 

solid angle 

Subscripts and superscripts: 

IN inclusive 

NN nucleon-nucleon 

P projectile 

T target 

Barred quantities represent matrices. 
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Figure 1. Inclusive proton scattering on 7 Li at 800 McV for scattering angle of 11° 



Figure 2. Inclusive proton scattering on 7 Li at 800 MeV for scattering angle of 15° 
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□ Experimental (ref. 16) 



Figure 7. Inclusive 4 IIe scattering on 1 He at 1.05 GeV/amu for scattering angle of 2.112°. 



Figure 8. Inclusive 4 He scattering on 'lie at 1.05 GeV/amu for scattering angle of 3.094°. 
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da/dp do, mb/(GeV/c)sr 



Figure 9. Inclusive 4 He scattering on 4 He at 1.05 GeV/amu for scattering angle of 3.630°. 



Figure 10. Inclusive 4 He scattering on 4 He at 1.05 GeV/amu for scattering angle of 4.552°. 
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Figure 13. Inclusive 4 He scattering on 16 0 at 1 GeV/amu for scattering angle of 4°. 



Figure 14. Inclusive 16 0 scattering on 16 0 at 1 GeV/amu for scattering angle of 0.5°. 
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Figure 15. Inclusive 1G 0 scattering on 16 0 at 1 GeV/amu for scattering angle of 1°. 


20 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503 

1. AGENCY USE ONLY(7eave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

May 1992 Technical Memorandum 

4. TITLE AND SUBTITLE 

Quasi-Elastic Nuclear Scattering at High Energies 

5. FUNDING NUMBERS 

WU 593-42-21-01 

6. AUTHOR(S) 

Francis A. Cucinotta, Lawrence W. Townsend, and John W. Wilson 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center 
Hampton, VA 23665-5225 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L- 17036 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM-4362 

I 11. SUPPLEMENTARY NOTES 1 


12a. DISTRIBUTION/AVAI LABILITY STATEMENT 

12b. DISTRIBUTION CODE 

Unclassified-Unlimited 



Subject Category 73 



13, ABSTRACT (Maximum 200 words) 

The quasi-elastic scattering of two nuclei is considered in the high-energy optical 
momentum transfer spectra for projectile ions are evaluated in terms of an inelastic 
corresponding to multiple knockout of target nucleons. The leading-order correction 
approximation is evaluated. Calculations are compared with experiments. 

model. Energy loss and 
multi pie- scattering series 
to the coherent projectile 






14. SUBJECT TERMS 

Nuclear reactions; Optical model; Inclusive reactions; Response functions 

15. NUMBER OF PAGES 

21 

16. PRICE CODE 

A 03 

17. SECURITY CLASSIFICATION 
OF REPORT 

18. SECURITY CLASSIFICATION 
OF THIS PAGE 

19. SECURITY CLASSIFICATION 
OF ABSTRACT 

20. LIMITATION 
OF ABSTRACT 


Unclassified 

roaw 


Unclassified 


Standard Form 298(Rev. 2-89) 

Prescribed by ANSI Std Z39-18 


298-102 


NAS A- Langley, 1992 























